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Abstract 

In this paper a Frenkel-Kontorova model with a nonlinear interaction potential 
is used to describe a vacancy defect in a crystal. According to recent numerical 
results [Cuevas et al. Phys. Lett. A 315, 364 (2003)] the vacancy can migrate when 
it interacts with a moving breather. We study more thoroughly the phenomenology 
caused by the interaction of moving breathers with a single vacancy and also with 
double vacancies. We show that vacancy mobility is strongly correlated with the 
existence and stability properties of stationary breathers centered at the particles 
adjacent to the vacancy, which we will now call vacancy breathers. 

Key words: Discrete breathers. Mobile breathers, Intrinsic localized modes. 

Vacancies, Breather-kink interaction. 

PACS: 63.20.Pw„ 63.20.Ry„ 63.50.+x„ 66.90.+r 



1 Introduction 



Discrete breatliers (DBs) are classical, spatially localized, time-periodic, nu- 
merically exact solutions which can be sustained by many non-linear lat- 
tices [1]. Their existence, which is proven by rigorous theorems [2] and a large 
amount of numerical results, is not restricted to special or integrable models. 
On the contrary, they can be found, in principle, in any discrete, nonlinear 
system and in any dimension [3,4]. They have been observed in experiments 
involving different systems, as Josephson-junctions arrays [5,6], waveguide ar- 
rays [7,8], molecular crystals [9] and antiferromagnetic systems [10]. They are 
also thought to play an important role in DNA denaturation [11]. 



Corresponding author. E-mail: jcuevas@us.es 



Preprint submitted to Physica D 



14 October 2005 




Fig. 1. Scheme of the equihbrium state of the Frenkel-Kontorova model with 
sine-Gordon substrate potential and Morse nearest neighbor interaction. The empty 
well represents a vacancy. 

Usually they are pinned to the lattice but under certain circumstances may 
become highly mobile [12]. In this case, an interesting problem arises: the in- 
teraction between moving discrete breathers (MBs) and local inhomogeneities. 
This problem has been addressed within different frameworks: scattering of 
breathers or solitons by impurities [13,14,15], by lattice junctions [16] and by 
bending points [17] of a given chain. 

Vacancies in a crystal are another type of local inhomogeneities. By using 
numerical methods, some of the authors found that a localized energy packet in 
the form of a moving discrete breather can put a vacancy into movement [18]. 
This vacancy migration induced by localized excitations had already been 
suggested as an explanation for some experimental results [19]. 

The aim of this paper is to find an explanation to the somehow qualitative 
results established in [18]. The main result is that vacancies mobility is highly 
dependent on the existence and stability of the breathers adjacent to the 
vacancy, hereafter called vacancy breathers (VBs). 



2 The model 



The simplest way to model a vacancy in a crystal consists in considering a 
one-dimensional chain of interacting particles submitted to a periodic sub- 
strate potential. It is known in other context as the Frenkel-Kontorova (FK) 
model [20]. A vacancy can be represented by an empty well of the substrate 
potential as Fig. 1 shows. 
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The Hamiltonian of this system is [15] 

N ^ 

H^Yl o^^l + + W{Xn - Xn-l) , (1) 

n=l ^ 



being Xn the absolute coordinate of the n-th particle. We have chosen the 
sine-Gordon potential 

V{xn) = ^[1 - cos(27ra:„/a)] (2) 



as the simplest periodic, substrate potential, with linear frequency normalized 

to the Tinity loq — ^JV"{0) = 1. We have chosen the Morse potential for the 
interaction between particles because its force weakens as the distance between 
particles grows and its hard part prevents the particles from crossing. It is 
given by 

W(Xn - Xn-l) = ^[e-'(-n-x„-^-a) _ i]2_ (3) 



The parameter C = W"{a) is the curvature of the Morse potential. We have 
taken C = 0.5 in order that MBs exist in this system for a breather frequency 
LUb — 0.9. The strength of the interaction potential can be modulated, without 
changing its curvature, by varying b, being a measure of the well width 
and C/2fe^ the well depth. We have also normalized the lattice period a and 
the masses to the unity. The dynamical equations are given by 

Xn + V'{Xn) + [W'{Xn - Xn-l) - W'{Xn+l - Xn)] = 0. (4) 



In this system, it is possible to generate DBs numerically using the standard 
methods from the anticontinuous limit [21]. We can also induce translational 
motion of DBs by using a simplified form of the marginal mode method [12]. 
It consists of adding a perturbation v — A(..., 0, — 1/-\/2, 0, l/\/2, 0, ...) to the 
velocities of the stationary breather, with the nonzero values at the neigh- 
boring sites of the initial breather center. The resulting DB kinetics is very 
smooth and resembles that of a classical free particle with constant velocity. 
Therefore, we can consider the total energy of a MB as the sum of a the inter- 
nal energy, equal to the one of the stationary breather, plus the translational 
energy, equal to the energy of the perturbation added K — 

In order to facilitate systematic studies, we have introduced dissipation at the 
boundaries. 
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(a) C=0.50, b=1, ^=0.1426 



(b) C=0.50, b=1, ^=0.1052 




35 r 

30 

25 



■^20 
O 



10 



20 



(c) C=0.50, b=0.8, ^=0.1367 
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(d) C=0.50, b=0.5, ^=0.1498 
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Fig. 2. Energy density plot for the interaction moving breather-vacancy. The particle 
to the right of the vacancy is located at n = 0. In (a), (c) and (d) the moving breather 
is reflected and the vacancy moves backwards, remains at rest and moves forwards 
respectively. In (b) the breather is transmitted and the vacancy moves backwards. 



3 Interaction of moving breathers with a single vacancy 



In order to investigate vacancy mobility, we have generated DBs far from a 
vacancy and then launched this breather against it. Our numerical calculations 
show that the outcome of the scattering is extremely sensitive to the initial 
conditions [18]. The incident breather can be reflected, trapped or transmitted 
(see Fig. 2), always losing energy as it occurs in the interaction between a 
MB and an impurity [13], whereas the vacancy can either move forward or 
backward or remain at rest. This scenario is very different to the one arising 
in the continuum limit where the vacancy (anti-kink) always moves backwards 
and the breather is always transmitted after the collision. 

We have studied the dependence of the vacancy mobility on the parameter 6, 
which controls the strength of the interaction between neighboring particles 
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Fig. 3. Probability for the vacancy to remain at its site (squares), move backward 
(circles) or forward (triangles). The interaction weakens when b increases. Note that 
no forward movement of the vacancy occurs to the right of the vertical line. 



of the chain. For each value of b we have analyzed a set of 601 numerical ex- 
periments of scattering corresponding to different initial conditions. As initial 
conditions we have considered a set of MBs with increasing kinetic energy, 
which were obtained choosing the values of the parameter A uniformly dis- 
tributed in the interval [0.10,0.16]. The lower bound of this distribution is 
chosen so that the threshold value for vacancies movement (see below) is sur- 
passed. The upper bound can be increased without any dramatic change in 
the results. However, an extremely high value of A can destroy the moving 
breather. 



Fig. 3 shows the probability that the vacancy moves forward, backward or 
remains at rest (averaged over the set of initial conditions) with respect to the 
parameter b. Note that, as it should be expected, the backward movement is 
the most probable behavior because in this case the incident breather pushes 
forward the particle to the left of the vacancy. In fact, the forward movement 
of the vacancy requires a strong enough interaction, i.e. b has to be smaller 
that a critical value bj ~ 0.7. 



Another interesting finding is shown in Fig. 4: the translational energy of 
the MB has to be higher than a minimal value K^in in a certain interval 
b G (0.5, 0.8) in order to move the vacancy. It is worth remarking that Kmin = 
means the minimum kinetic energy needed to move a breather, which is 
actually different to zero. However, we have included this notation in order to 
have a clear picture of the vacancy movement. 
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(a) Single vacancy 
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Fig. 4. Minimal translational energy (Kmin) to move a vacancy. 
4 Vacancy modes and vacancy breathers 

We can try to understand these resuhs by studying: a) the hnear modes, which 
can be obtained by linearizing Eqs. 4 around the equihbrium state shown 
in Fig. 1; and b) the vacancy breathers, i.e., the nonhnear locahzed modes 
centered at the particles adjacent to the vacancy, with the same frequency as 
the incident MB. 

Fig. 5a shows the dependence of the linear frequencies with respect to b. Note 
that the minimum linear frequency corresponds to a linear localized mode 
which we call linear vacancy mode (LVM). On the other hand, there are three 
types of VB: 1-site vacancy breathers (IVB), which consists of a single excited 
particle adjacent to the vacancy; and 2-site vacancy breathers, obtained by 
exciting the two particles adjacent to the vacancy, which can vibrate in-phase 
(2VBp) or anti-phase (2VBa). Fig. 5b shows the bifurcation pattern for these 
VBs. A letter at the beginning of the abbreviation specifies the breather sta- 
bility: stable, /i [/ harmonically unstable, sU subharmonically unstable or oU 
oscillatory unstable. For example sU2VBa means a subharmonically-unstable, 
two-site, vacancy breather in anti-phase. 

For weak interaction (6 large) there exist stable one-site VBs at both sides of 
the vacancy as stated by the theory [1]. As b decreases the SIVB disappears 
through an inverse pitchfork bifurcation to two unstable two-site VB in phase 
(hU2VBp) which becomes stable (S2VBp). Further decrease of b brings about 
the annihilation of the S2VBp when the breather frequency coincides with the 
frequency of the LVM with the same profile. 

The two-site VB in antiphase is stable (S2VBa) for the largest values of b 
represented in Fig. 5b. When the interaction increases [b decreases) this 2VB 
becomes oscillatory unstable (oU2VBa), and finally subharmonically unstable 
(sU2VBa). 
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Prom this bifurcation diagram we can deduce that the stabihty of the VBs 
is correlated to the existence of a kinetic energy threshold to move the va- 
cancy, in the interval b e (0.5, 0.8). Comparison of Fig. 4 and 5b suggests that 
the existence of a subharmonically or harmonically unstable, 2- site, vacancy 
breather is a necessary condition in order to get an optimum vacancy mobility 
(no kinetic energy threshold). 



Generally speaking, any approaching breather with any kinetic energy is able 
to move the vacancy if at least one the two 2VBs exists and is either sub- 
harmonically or harmonically unstable. We cannot expect, however, an exact 
agreement between the numerically exact bifurcation values for a given fre- 
quency with the observed changes in behavior for the simulations. The reason 
for that is that a moving breather is obtained by perturbation with an asym- 
metric mode. Therefore, it is no longer an exact solution of the dynamical 
equations, its frequency is shifted and it is not unique. The bifurcation values 
for other frequencies will change and interaction of phonons emitted by the 
MB are expected to have an (unknown) influence, hence the probability anal- 
ysis of the simulations outcome. However, the only important exception to the 
first assessment in this paragraph in the region b G (0.73,0.8) is noteworthy. 
The role of the unstable vacancy breather is probably to act as an interme- 
diate structure with high oscillations for the particles nearest to the vacancy, 
which, being unstable, leads to one of them changing place. Clearly the stable 
S2VB does not play that role and its existence in the vicinity of that region 
seems to be the reason for that discrepancy -proximity in b for a = 0.9 also 
means proximity in a;b for the a given value of b, which can be excited by the 
bundle of the MB frequencies-. 



Therefore our conclusion is that vacancy mobility is governed, for a particular 
choice of the substrate potential, by the strength of the interaction between 
the particles adjacent to the vacancy, which is correlated to the existence and 
stability properties of vacancy breathers. This correlation manifests in two 
different ways. Firstly, the interaction has to be strong enough to move the 
vacancy forward, or equivalently, there must exist no linear vacancy modes. 
Secondly, a threshold value of the translational energy of the incident MB in 
order to move the vacancy does not exist if an harmonically or subharmonically 
unstable VB exists. 



Another choice of the interaction potential would alter the values of b which 
separate the different regimes. However, the correlation between the vacancy 
mobility and the strength of the potential would be similar. 
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Single vacancy. C=0.5 




Fig. 5. a) Linear spectrum. Also represented is the breather frequency = 0.9. 
Note the existence of a locahzed Hnear mode, (b) Bifurcation diagram for the non- 
Unear locahzed modes obtained by exciting the particles adjacent to the vacancy. 
The bifurcation variable A is the difference between the relative displacements of 
the particles to the left and to the right of the vacancy. The vertical lines indicate 
bifurcation points. See the text for the breather codes. 

5 Double vacancy 



We have studied the mobihty of a double vacancy (two empty neighboring 
wells of the substrate potential) to test our conjecture. The outcome is similar 
to the single vacancy case except for the fact that the left vacancy never moves 
forwards. Thus, if the right vacancy moves forwards, the double vacancy splits 
in two (see Fig. 6). Fig. 7 shows the linear spectrum and the bifurcation 
diagram for the VB in this case. In these figures the lowest value of b is 
approximately 0.5 since there are no equilibrium states for b < 0.5. This 
result also implies the non existence of double vacancies equilibrium states 
for a harmonic interaction potential. As b increases, we observe two-site VBs 
emerging from linear localized modes, and afterwards an inverse pitchfork 
bifurcation similar to the one described for a single vacancy. There is not 
any harmonically or subharmonically unstable VB for b < 0.75. Therefore, 
we should expect that there is a minimum value for the kinetic energy of the 
incident MBs able to move the double vacancy. This is confirmed by numerical 
computations as Fig. 8 shows. 



6 Summary 



We have observed numerically the interaction of moving breathers with va- 
cancies in the simplest, physically consistent model. We have found that the 
breathers can be reflected, transmitted and trapped, both by the single and 
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Fig. 6. Energy density plot for the interaction moving breather-double vacancy, in 
the vacancy splitting regime. The particle to the right of the vacancies is located at 
n = 0. 



Double vacancy. C=0.5 




Fig. 7. Double vacancy case: a) Linear spectrum. The breather frequency Wb = 0.9 
is also represented. Note the existence of two localized linear modes, (b) Bifurcation 
diagram for the nonlinear localized modes obtained exciting the particles adjacent to 
the double vacancy. The bifurcation variable A is the difference between the relative 
displacements of the particles to the left and to the right of the double vacancy. The 
vertical lines indicate bifurcation points. See text for the breather codes. 

the double vacancy. The single vacancy can move forward, backwards and re- 
main and rest, but the double vacancy cannot move forward, instead it can 
split into two single vacancies, one at rest and the other moving forwards. 
This phenomenology is very different from the continuous analogue, where 
the breather is always transmitted and the vacancy (anti-kink) always moves 
backwards. 

The vacancy migration is strongly correlated with the existence and stability 
of vacancy breathers, which we have studied in detail. In order to move the 
vacancy the interaction has to be strong enough and the energy of the incident 
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(b) Double vacancy 
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Fig. 8. Minimal translational energy (^K-min^ to move a double vacancy. 

MB has to be above a threshold. This threshold disappears approximately 
around bifurcation points where vacancy breathers become unstable. 

We think that the study of the properties of breathers next to a defect may help 
to understand the mobility of other class of point defects such as interstitials, 
which might be related with recent defect migration observed experimentally 
in ion-irradiated silicon. 



7 Acknowledgments 

The authors acknowledges sponsorship by the Ministerio de Educacion y Cien- 
cia, Spain, project F1S2004-01183. They also acknowledge JC Eilbeck and F 
Palmero for useful discussions. 



References 

[1] S Flach and CR Willis. Discrete breathers. Phys. Rep., 295:181, 1998. 

[2] RS MacKay and S Aubry. Proof of existence of breathers for time-reversible or 
Hamiltonian networks of weakly coupled oscillators. Nonlinearity, 7:1623, 1994. 

[3] Localization in nonlinear lattices. Physica D, 119:1-238, 1999. Focus issue 
edited by S Flach and RS Mackay. 

[4] Nonlinear localized modes: physics and applications. Chaos, 13:586-799, 2003. 
Focus issue edited by Yu S Kivshar and S Flach. 

[5] P Binder, D Abraimov, AV Ustinov, S Flach, and Y Zolotaryuk. Observation 
of breathers in Josephson ladders. Phys. Rev. Lett., 84(745), 2000. 



10 



[6] E Trias, JJ Mazo, and TP Orlando. Discrete breathers in nonlinear lattices: 
experimental detection in a Josephson array. Phys. Rev. Lett., 84:741, 2000. 

[7] HS Eisenberg, Y Silberberg, R Morandotti, AR Boyd, and JS Aitchison. 
Discrete spatial optical solitons in waveguide arrays. Phys. Rev. Lett., 81:3383, 
1998. 

[8] JW Fleischer, T Carmon, M Segev, NK Efremidis, and DN Christodoulides. 
Observation of discrete solitons in optically induced real time waveguide arrays. 
Phys. Rev. Lett, 90:023902, 2003. 

[9] BI Swanson, JA Brozik, SP Love, GF Strouse, AR Bishop, W-Z Wang, and 
MI Salkola. Observation of intrinsically localized modes in a discrete low- 
dimensional material. Phys. Rev. Lett., 84:741, 2000. 

[10] UT Schwarz, LQ English, and AJ Sievers. Experimental generation and 
observation of intrinsic localized spin wave modes in an antiferromagnet. Phys. 
Rev. Lett, 83:223, 1999. 

[11] M Pcyrard. Nonlinear dynamics and statistical physics of DNA. Nonlinearity, 
17:R1-R40, 2004. 

[12] Ding Chen, S Aubry, and GP Tsironis. Breather mobility in discrete lattices. 
Phys. Rev. Lett, 77:4776, 1996. 

[13] J Cuevas, F Palmero, JFR Archilla, and FR Romero. Moving discrete breathers 
in a Klein-Gordon chain with an impurity. J. Phys. A: Math, and Gen., 
35:10519, 2002. 

[14] K Forinash, M Peyrard, and BA Malomed. Interaction of discrete breathers 
with impurity modes. Phys. Rev. E, 49:3400, 1994. 

[15] OM Braun and YuS Kivshar. Nonlinear dynamics of the Frenkel-Kontorova 
model. Phys. Rep., 306:1, 1998. 

[16] I Bcna, A Saxcna, and JM Sancho. Interaction of a discrete breather with a 
lattice junction. Phys. Rev. E, 65:036617, 2002. 

[17] GP Tsironis, JM Sancho, and M Ibahcs. Localized energy transport in 
biopolymcr models with rigidity. Europhysics Letters, 57:697, 2002. 

[18] J Cuevas, C Katerji, JFR Archilla, JC Eilbeck, and FM Russell. Influence of 
moving breathers on vacancies migration. Phys. Lett. A, 315:364-371, 2003. 

[19] P Sen, J Akhtar, and FM Russell. Mev ion-induced movement of lattice disorder 
in single crystalline silicon. Europhys. Lett., 51:401, 2000. 

[20] Yal Frenkel and T Kontorova. Phys. Z. Sowietunion, 13:1, 1938. 

[21] JL Marm and S Aubry. Breathers in nonlinear lattices: Numerical calculation 
from the anticontinuous limit. Nonlinearity, 9:1501, 1996. 



11 



